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Abstract 

In this work we study possibility of chaos formation in the dynamics governed by paradigmatic 
model of Cavity Quantum Electrodynamics, the so called James-Cammings model. In particular 
we consider generalized JC model. It is shown that even in the case of zero detuning dynamics is 
chaotic. Kinetic approach for the problem under study has been applied. 



PACS numbers: 73.23.-b,78.67.-n,72.15.Lh,42.65.Re 



FIG. 1: Schematic representation of atomic levels and available inter level transitions for generalized 
JC model. 



INTRODUCTION 

Cavity quantum electrodynamics (CQED) is a rapidly developing field ofphysics studying 



the interaction of atoms with photons in the high-finesse cavities 4j. Interest to 

such a systems basically is caused by two facts: One of them is the possibility of more deep 
understanding of quantum dynamics of open systems. Second argument is the possibility of 
practical application in the field of quantum computing [^j]. 

In particular CQED experiments implement a situation so simple that their results are of 
great importance for better understanding of fundamental postulates of quantum theory {(| . 
They are thus appropriate for tests of basic quantum properties: quantum superposition, 
complementarily or entanglement. In the context of quantum information processing, the 
atom and cavity are long-lived qubits, and their mutual interaction provides a controllable 
entanglement mechanism an essential requirement for quantum computing Q|. 

In general dissipation processes must be taken into account when discussing problems of 
CQED. In particular there are two dissipative channels for system: the atom may sponta- 
neously emit into modes other then preferred cavity mode, and photons may pass through 
the cavity output coupling mirror. 

But modern experiments in CQED have achieved strong atom-field coupling for the 
strength of the coupling exceeds both decay processes si, 9, 10]. If so, then problem is 
reduced to the driven Jaynes-Cummings (JC) Hamiltonian, which models the interaction of 
a single mode of an optical cavity having resonant frequency with a two level atom comprised 



2 



of a ground and exited states 11]. 

This model is basic model of interaction of radiation with matter and describes the energy 
exchange between atom and quantized radiation field in an ideal lossless cavity. 

In further was offered JC model generalized for three level optical atoms (see FigHT). For 

nn 

more details see ill njl. 

In most general case atom-radiation field interaction should involve not only the internal 
atomic transitions and field states but also the center-of-mass motion of the atom. With the 
recoil effect taken into account, generalized JC Hamiltonian takes form 

2 3 2 
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H = — + fiu<Mb a + S ^E j R jj + h^2g a (b a R 3a + b+R a3 ) cos(a%J (1) 

a=l j=l o=l 

here x and p are coordinate and pulse of the atom R 3a , R a3 are transition operators 

3 

= 1,2 , Rjj is the operator of the level population, satisfying the condition ^2 Rjj = 1. 



In the basis of atomic states 
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Rjj operators are the generators of the SU (3) 
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Model given by the Hamiltonian (1) is nonlinear and non-integrable and due to this 
may display chaotic behavior. Analogous problem but for ordinary two level JC model 
was studied in [jjj] in detail. Result obtained in is that dynamics essentially depends 
on detuning parameter between inter level transition frequency and frequency of radiation 
field 5 = u a — E jfi. Namely, in case of zero detuning system displays regular behavior. 
Goal of this work is the study of nonlinear dynamics governed by Hamiltonian (1). For 
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better understanding of the problem we shall discuss resonant case. We shall try to find out 



difference in dynamics between generalized and ordinary JC model studied in 

In resonant case, complete set of Heisenberg equations of motion corresponding to the 
Hamiltonian (1) looks like 

T = otp(r), 

^f-^ = u\(t) sin(x) + ft ku 2 (r) sin (for), 
= z/ i( r ) cos ( x ), 

dR 2 (T) 



ft v 2 {j) cos(kx), 
- 2cos(x)(M 1 + 1)(1 - 2Rx(r) - R 2 (t)) - ncos(kx) ■ B(t), 



(It 
dr 



(2) 



= 2ftcos(A;x)(M 2 + 1)(1 - R^t) - 2R 2 (t)) - cos(x) • B{r) 
^ = ncos(kx)-C(r), 
*gd = -oos(fcr).C7(r) ) 

^ = cos(A;x)(M 1 + l)z/ 2 (r) + ftcos(£;x)(M 2 + l)i^(r), 
^ = cos(x)(Mi + 1)u 2 (t) - ftcos(A;x)(M 2 + l)u x {r) 
—fir 1 = cos(x) • i/i(r), 

= ft cos(fcr) ■ i/ 2 (r) 

In (2) transform to the Hermit variables A\ = i(biR 3 i — bfRis), A 2 = i{b 2 Rz2 — b 2 R 2 z), 
Ui = (kR 3 i + btRu), U 2 = (b 2 R 32 + btR 23 ), B = b x biRn + bfb 2 R 12 , C = i(bib+R 21 - 
bfhRu), Ni = bth,N 2 = b+k, 



and procedure of semi-classical averaging is done jl^] x = fc/j < x >, p 



<p> 

hkt 



k f2 



-Ri =< #n >, R2 =< R22 >, i>i =< Ax >, =< A 2 >, Ul = (u^, u 2 = (u 2 y k - 
B = C = (c^, t = g±t, a = ft = ^ . In addition following notations for 

motion integrals are introduced Mi = Ni — Ru = const, M 2 = N 2 — R 22 = const. 

Below the set of equations (2) will be the object of our interest. We shall consider two 
cases: 

1) Variable x is slow as against Ui(r), u 2 {t). In this case set of equations (2) separates 
into two subsystems 
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^ = ncos(kx) -C(r) 
*^ = -oo8(x).C(r) l 



= cosOXMa + l)u 2 (r) - ficos(A;x)(M 2 + l)«i(r), (3) 

c±e(t) 



- iii(r) sin(x) + Vtku 2 {r) sin(fcx). 



and 



= ^(r) cos(x), 
^=n^(T)cos(fcr), 

- 2cos(x)(M 1 + 1)(1 - 2R l {r) - R 2 (r)) - ft cos(kx)B(r) 



dr 

dr 
dB{r) _ 

dr 
djVi(r) 

dr 



2ttcos(kx)(M 2 + 1)(1 - Rx(t) - 2R 2 (r)) - cos(kx)B(r), (4) 
cos (x) (Mi + l)r/ 2 (r) + Qcos(kx)(M 2 + l)^i(r), 
= cos(x) ■ ^i(r), 



^rl = ncos(kx)-u 2 (r), 
After solving of first three equations of (3), we have 



C(t) = Aan(n,r), 

Ul (r) =u 1 (0)-^cos(n 1 r), (5) 



M2 ( r ) = M2 (0) + ^cos(fi 1 r) 



where 



^2 = C0S 2 (x)(M 1 + 1) + fi 2 cos 2 (fcr)(M 2 + 1). 
Then if, k = ^ = 1 , motion of the atom inside of quantum cavity satisfies the equation 



A' 



d 2 x 



2 + |a|(ui(0)+« 2 (0))sin(z) =0, (6) 



with the possible solutions 



Jcra(r,e), 0<1 



p = 2Qu,l , (7) 

[dn(r, 1/0), 0^1 



where = ± (1 +H/w%), = |a|(tti(0)+tt 2 (0)), ^I+u; 2 , cos(x ) and cn(...), 

are Jacobi elliptic functions 15] . Two solutions (7) correspond to the two different phase 
trajectories. One of those solutions dn(...) corresponds to the closed phase trajectory. Other 
one corresponds to open phase trajectory. In case of closed phase trajectory it means that 
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after the time interval , atom iterates initial state. For solving of the subsystem (4) we can 
use the method given in [yj]. As a result we get: 



Ri(t) = /i(cos(Ar) - 1) +/?sin(Ar) + A 2 [u(cos(2Ar) - 1) + vsin(2Ar)] + #i(0), 
R 2 (r) = -/i(cos(Ar) - 1) - /3 sin(Ar) + A|[u(cos(2At) - 1) + v sin(2Ar)] + R 2 (0) 

where 



(8) 



ix = X-^XiRxiO) - XiR 2 (0)] + (Xi - Xf)K], A = y Af + A|, 
= A" 3 [A'i?i(0) - A 2 i? 2 (0)], u = ^A~ 4 [A 2 [2i?i(0) + 2i? 2 (0) - 1] + K], 
u = Ix^lR^O) + R 2 {0)}, Ai = cos(x) v / Mi + l, A 2 = Q cos(kx) ^ M 2 + 1, 



and 



N±(r) = /i(cos(Ar) - 1) + /5sin(Ar) - A?[w(cos(2Ar) - 1) + z/sin(2Ar)] + 7Vi(0), 

N 2 (t) = -^(cos(Ar) - 1) -/3sin(Ar) + A^[n(cos(2Ar) - 1) + v sin(2Ar)] + N 2 (0). 

Really the quantities Ri(r), R 2 (t), Rs{t) = 1 — R\{r) — R 2 {r) are observable on the 
experiment. But in general case, when variable x(r) is not adiabatic, we have to investigate 
(2) by use of numerical methods. If x(r) is non-adiabatic variable, origination of chaos in 
the system is quite possible. In this case dynamical consideration loses sense and transition 
from dynamic to the statistical description is needed. As the result, statistical conceptions, 
namely Kolmogorov entropy and fractional dimension becomes important. Those concepts 



being statistical are strongly related to t 



local instability of phase trajectories 



re systems motional characteristics, namely to the 
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Results of numeric calculations, for the 



realistic values of parameters corresponding to the high-finesse Fabry-Perot cavity and real 
atoms 0, Q are presented on Figj2j FigJHl 

For determination of the width of Furrier transform of the correlation function we used 
method of fast Furrier transform 

Gxij) =< x(t + r)x(t) > 



G x (u) = I drG x (T)exp(tujT) 



1 + w 2 r c 2 
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FIG. 2: Plot of atom mass center coordinate dependence on time. As one can see from this Fig. 
motion is chaotic. As every other numeric result, plot is obtained for the values of parameters 
a = 0.01, = 0.5, k=3, M 1 = 1.6, M 2 = 1.7 



T 

where < (...) >= lim i J\...)dt means time average, r c is the correlation time. Result 

T->oo 

of numeric calculations are presented on FigJHl Nonzero width of the Fourier form is the 
sing of chaos. 

One more sign for confirmation of chaos existence is the values of maximal Lyapunov 
exponent. In this sense it is worth to compare our result with the result obtained for 
ordinary JC model [lJ|. For numerical calculation of maximal Lyapunov exponent we shall 



use G. Benettin's algorithm consisting in the following procedures [18[]: Solving two different 
systems starting from the initial points xq, xq + x®, \x~q\ = e, where e is small initial distance 
between phase trajectories. Then after each step x(ti) = Xi, x(ti) = X{ one has to do re- 
scaling x\ = exi/\xi\ use as initial conditions x i; Xi + and so on up to the xn,%n + ■ 
Then Lyapunov exponent is given via 



A 



1 N 
N ^ 

k=l 



e 



(9) 



Result of numeric calculations for the maximal Lyapunov exponent is presented on Fig.4. 



It is evident from the FigJH that in case of ordinary JC model corresponding to the values 
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FIG. 3: Fourier transform of the correlation function 

oo 

p(to) = G p (u>) = J dr exp(zo;r)Gp(r). Nonzero width of the correlation function is the sign of 
o 

classical dynamical stochastisity 




FIG. 4: Maximal Lyapunov exponent as a function of ratio between constants of interaction between 
atom and radiation field (f2 = 52/51) 



of parameter Q = 0, maximal Lyapunov exponent is equal to zero A = 0. This result is in a 
good agreement with the results obtained in 

Other important physical quantity characterizing classical dynamic stochasticity is the 



S 



fractional dimension of the systems phase space. 
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or defining of fractional dimension we 
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2l|. 



will use algorithm of P.Grassberger, I.Procaccia 

Let we have set of state vectors (xi, i = 1, 2, • • • , N) corresponding to the successive steps 
of numeric integration. In our case x% is t ne complete set of variables (2) with the values 
corresponding to the moments of time t = tj. Then we can use numeric data for estimation 
of following expression 



C(e) = lim — — — — 



N 



\Xi-Xj\) 



where 9{x) 
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x < 

1 x > 



is the step function. According to the P.Grassberger, I.Procaccia 



2l| fractional dimension may be defined as 

C(e) 



D = lim 



e -o bg(e) 

Expected dependence of C(e) is e D . So, plot must be a line with angular coefficient D. 
Result of numeric calculations are presented on Figj5j 



II. QUANTUM MECHANICAL CONSIDERATION. MIXED STATE FORMA- 
TION. 

On the basis of numeric results of previous section, we can conclude that dynamics is 
chaotic. In case of quantum consideration we neglect kinetic energy of atomic motion as 
compared to the level transition frequencies, and random nature of motion we shall try to 
take into account by considering x(t) as time dependent random Wiener process. As a result 
in interaction representation we get: 

M(t)) 



dt 

where interaction operator is 



v\m) (io) 



V = g a (b a R 3a + b+R a3 ) ■ cos k fa x (11) 

a=l 

Let present initial wave function as a direct product of atomic and field states 

1-0(0)) = l^atom) ® \lp field) 




0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 



FIG. 5: Dependence of C(e) on the values of e , plotted by numerical integration of the set of 
equations (2) , for the following values of parameters a = 0.01, Q, = 0.5, k = 3, M\ = 1.6, M2 = 1.7. 
A solid line corresponds to least-squares approximation of the results of data processing. According 
to this plot, fractal dimension of the system D = ^^(gaj-inCei) 1 ^ * s ec l ua l to Z) ~ 2.44 



where 



\lp field) = W m\ n l) + Yl W ^\ n 2), \lf> atom) — Cai| a l) "I" C a2 \a,2) + C a3 \a^) 

ni=0 712=0 



ai 








\a 2 ) 




1 

V / 



\ a s) 





V / 



As far as operator (11) can mix only following states |a 3 ni), \a± ni + 1), |a3 77-2), 1 0-2 ^-2 + 1)> 
solution of (10) we will search in the form 



(12) 



CLjTli 



After substitution (12) in (10) and solving equations for C a . ni {t) we get 



Casm (t) n 1 6 



Ci, 
2 



J coskf x(t')dt' iXi f cos kf x(t')dt' 

+ ^fe ° 



„ — iAi f cos kr, x(t')dt' _ iAi f cos k t, x It' )dt' 

C ..(t) = 1 _ Q2. e 1 

w ami+l Vv 2 2 



(13) 
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where Ai = giy/rfi + 1 



C\ — C a . jni (0) + C aini +i(0); C*2 — C aani (0) — C aini+ i (0); 

q —i\2 /cos kf 2 x(t')dt' q i\2 J cos kf 2 x(t')dt' 



a 



(t) 



—i\2 J coskf x(t')dt' 1X2 J coskf x(t')dt' 

±dLp — k±p 



(14) 



Here A 2 = 92V n 2 + 1 



C3 — C a3Tl2 (0) + C a2n2+ i(0); C4 — C a3ri2 (0) — C a2n2+ i(0); 



Density matrix of the system: field + atom is direct product of two density matrixes 



Pijkl — Patom ® P field ~ C ainj • C* 



a k n t 



(15) 



where (...) means averaging. 

Sign of mixed state formation is zeroing of non-diagonal matrix elements of density matrix 
(15), 22, l23|. This process is non-reversible because, along with zeroing of non-diagonal 
matrix elements, loss of information about wave function's phase factor takes place. So, 
to prove irreversibility of dynamics, it is sufficient to show zeroing of non-diagonal matrix 
elements. 

Substituting (13), (14) into (15), it is easy to see that non-diagonal matrix elements 
contain terms like 



C ai nMC* (0)exp 



z J u)f 2 (t')dt' 
. 



(16) 



where tof 2 (t') = g\ cos(fc/ 1 x(t')) ± #2 cos(kf 2 x(t')). Due to the random nature of x(t'), below 
we will consider ujf 2 as a random process and for simplicity omit indices. It is clear that 
exponent in (16) may be considered as a functional of the random function u(t) 



Q[u] = exp 



i J u(t')dt' 
. 



(17) 



So, to get non-diagonal matrix elements of the density matrix (15), we need to do statistical 
averaging with respect to the all possible realizations of random function u(t) 



p.M!%~ c ^(o)c; ni (°) (QMt]]) 



(18) 
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where ((. . .)) means statistical average. Average values of the functional (18) may be calcu- 
lated by doing following continual integral 



Q[uj] = exp 



i j (jj(t')dt' 
L o 



lim / dujN ■ ■ ■ doji exp 



k=l 



Pn[u] (19) 



where At k = t^ k ' — t^ k x \ = 0, t^ N ' = 1, -P/vM is the multi-dimensional normal distribu- 
tion function given by 



P N [u] = (2n)' N J d\ 1 ...d\ N e*v 



d\x exp 


-i^X k u k 


exp 




k 





9 ^2 C kk'^k\' 



(20) 



where X k are the distribution parameters and C kk > is the covariation matrix [24] 
After substituting (20) in (19) and doing integration we obtain 



N 



/ duJi . . . duj N exp[i ^ UkAt^P, 



k=l 



f d\\ . . . (iAjv exp 



k,k' 



N 

x II / exp[iu k (At k - X k )}du k } 
k=i 



(21) 



= fd\ 1 ... dX N 5(X 1 - At 1 )6(X 2 - At 2 ) ...5{X N - At N ) x exp 
So, for average value of the functional we get 

i£c(t«^'))At fc At fe , 



\ £ At k At k , 



{Q[oj}) = lim exp 



k,k 



exp 



k,k' 



t t 



- / dt j c(t,t 







(22) 



For normal random process c(t , t") = c(t — t ). Then making transform to the new variables 
t' — t" = r, t' + 1" = £ and doing integration in (22) over variable £ we get 



(Q[u]) = exp 



t / drc(r) 



(23) 



After assuming that correlation function c(r) = (u;(t + r)u)(t)) has the Gaussian form c(r) 



from (23) we get 



(Q[u\) « exp 



- [Z* Erf[t^\ 



(24) 



where Erf(. . .) is the error function 24] 
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So, for non-diagonal matrix elements we have 

irt 



■Erf[t^\ 
a 2 



(25) 



From (25) it is clear that after the lapse of time t > ^/oqJtx zeroing of non-diagonal 
matrix elements and formation of mixed state happens. 

III. KINETIC DESCRIPTION 

In previous section we investigated mechanism of mixed state formation. After mixed 
state is formed, quantum-mechanical consideration loses meaning and there is a need to use 
kinetic description. Kinetic equation for chaotic quantum-mechanical system first time was 
obtained in 16|. But this study was done in the semi-classical domain. Namely, in zero- 
ing of non-diagonal part of density matrix was proved by use of semi-classical approximation. 
Our purpose is to do the same in exceptionally quantum domain, without application to the 
semi-classical methods. As the zeroing of non-diagonal part is already shown, problem is 
partly solved. 

It's clear that system atom+ field is adiabatically isolated. So entropy product is zero. 
But by considering field as a thermostat, we can calculate entropy growth for atomic sub- 
system. If AS > 0, then process in non-reversible. 

As we have already mentioned at t — 0, atom is not connected with light and state vector 
\ip{t = 0)) is the direct product of two states 

foKO)) = |VWn> ® Afield) (26) 

where 

oo oo 

l^fieid) = ^2 Wm\ni) + ^2 W n2\ n 2), \ipatom) = C ai \a x ) + C a2 \a 2 ) + C a3 |a 3 ) (27) 

n\=0 U2=0 

Let us note that interaction (11) mixes only states 

\a 3 m), \ai m + 1), \a 3 n 2 ), \a 2 n 2 + 1). 
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Then comparing (26) with (12) one can see that initial conditions in (13), (14) are of form 

Ca 3 m(0) = W ni C aa , 
Cami+i(0) = W ni+1 C ai , 

Ca 3 n 2 (0) = W n2 C a3 , 
Ca 2 n 2 +l(0) = W n2+ iC a2 . 

Let us remind that inversion of atomic populations / is a value of interest, since it can be 
measured on the experiment. Inversion is defined as difference between level populations 



L a s a 2 



(t) = W(t,\a 3 ))-W{t,\a 2 )), 

Iasatit) = W(t,\a 3 ))-W(t,\ ai )), 



where 



W(t,\a))=Y J \Can(t)\ 2 . 
n=0 



(29) 



(30) 



After substituting (13), (14) in (30), and using initial conditions (28) we get 



W(t, |a 3 » = \ E \(W n2 C as + W n2+1 C a2 )Q* M + (W n2 C a3 - W n2+1 C a2 )Q n2 [uj 2 }\ + 



rt2=0 



+i £ \(w ni c a3 + ^ 1+1 c ai )Q; l H + (w ni c a3 -w ril+1 c ai )Q ni H| 2 , 

m=o 

oo „ 

W(t, |a 2 )) = \ E \{W n2 C a , + W n2+1 C a2 )Q* n2 [u 2 ] - (W n2 C a:j - W n2+1 C a2 )Q n2 [u 2 }\ 2 , 

rt2=0 

oo „ 

w(t,\ ai )) = \ E |(^n 1 a 3 + ^„ 1+1 c , a jg;j^ 1 ]-(^ ni c , a 3-w ni+1 c , a jQ n j^ 1 ]|\ 



ni=0 



To get we also have introduced following notation for unitary functional 



Qm, 2 [^1,2] = exp 



-%,2V n i,2 + l y cos(k fl 2 x(t))dt 





(31) 



(32) 



q; 1i2 K 2 ] = Q~l 2 M, Q ni MdQnlM = I- (33) 

By using of (31) it is easy to define inversion values for any moment of time. But before 
doing this let us mention one interesting fact. In ordinary JC model which for zero detuning 
is not characterized by chaos, inversion is characterized by periodical revivals in time. See 
for example 25[ . So absence of periodical revivals maybe useful for experimental observation 
of quantum chaos. After using (32), (33) from (31) we get 

oo 



n 2 =0 



(34) 
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(W(t,\ai))) = E ^(<^V<+i^)-j(<^V< + i^ 2 J((^[^]> + (^ 1 2 [^ 

ni=0 

(35) 

where (...) again means statistical average. Then taking into account (24) 
from (34), (35) we get 

OO 

(W(t, |a 2 »> = £ [\(Wt 2 Cl + Wl 2+1 Cl 2 )l 

U2=0 

oo 

l«i»> = E + W^ 1+1 C^)], (36) 

711=0 

(W{t,\a 3 ))) = (W(t,\ ai ))) + (W{t,\a 2 ))) . 
At last let us remind that W% describes field states and obeys to Poisson distribution 



26] 



-"1,2 / - \ 

„ n,, expl — n-\ 2 ) 

and C 01 , C a2 , C az are probabilities of level occupations. 

Let us assume that, at the initial moment of time, atom is in the lowest state 

p ai (< = 0) = |C a mi+i(t=o)| 2 = 1, Pa2 (t = 0) = Pa3 (t = 0) = 0. (38) 

Then initial values of entropy is zero 

3 

S(t =0) = E PaAt = 0) HPaAt = 0)) = 0, (39) 

di=l 

so system is in pure quantum-mechanical state. 

After the lapse of time t = to more than time of inter level transition to ~ l/fl'aj system 
may perform multiple transitions between levels. That is why probability to find system in 
other states will be nonzero: 

C ai ^0, C a2 ^0, C a3 ^0, t>t . (40) 

Despite of this fact to talk about probability of population of different states is early yet. 
The point is that in time interval: 

to < t < (41) 

V 71 

interferentional terms are nonzero. Therefore the state of the system will be pure one. But 
unlike of the initial state (38), which is simple state, the state of the system in time interval 
(41) is superposition one. 

15 



Superposition state is pure quantum mechanical state and only after zeroing of interfer- 
entional terms in (18) superposition state passes to mixed one. Such a transition occurs in 
times: 

t > M (42) 

But in time interval (41) while the system is in pure superposition state, from the symmetry 
point of view, it is clear that the coefficient values(40) have to satisfy the following relation: 

C ai (t <t< yj^j ~ C a2 (t <t< ~ C a:i (t <t< ~ C. (43) 

Before proceed to the more exact determination of the quantities (36), we note one important 
fact. Values of the parameters (28) C ai , C a2 , C a3 , taken at the initial moment of time t — 0, 
define level populations. At the moment of time t = atom and cavity field are not connected 
with each other (26). Therefore the following relation should be hold: 

C 2 ai +Cl 2 +Cl = l. (44) 

After the laps of time interaction (11) will mix atom and field's states and due to the fact 
that (41) does not correspond to the normalization condition of complete wave function, it 
will not be valid henceforth. Considering of level populations again is possible only after the 
time interval (42). But in this case level populations are defined by the quantities (36), (28) 
and not by the coefficients (69). fficients (69). For determination of level's populations we 
shall use (36), (43) and take into account that fields states W ni2 obeys Poisson distribution 
(37). As a result we obtain 

Then taking into account (36) result we get 

112 

All = 4, Pa 2 = 4, Pa 3 = -, (45) 

Result is very interesting. Most populated is high exited level p a3 . This means that we have 
stochastic absorption of field energy. For entropy growth we have 

AS = -(In4 + ln2) > (46) 

So, process is non-reversible. 

IV. CONCLUSIONS 

Let sum up and analyze the results obtained in conclusion. 
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The aim of this work was to study generalized JC model being subject to resonator field. 
Interest to such a systems is caused by the fact that they are the most perspective to be 
used in quantum computer. The question that came up is the following: by how much will 
be state of the system controllable and dynamics reversible? We have considered the most 
general case, when interaction of the system with field depends on coordinate of the system 
inside resonator. 

Contrary to generally accepted opinion, it has turned out that the absence of detuning, 
between resonator field and frequency of the inter-level transitions, does not guarantee re- 
versibility of the system's state. During evolution in time the system executes irreversible 
transition from pure quantum-mechanical state to mixed one. At the same time, the time 
needed for formation of mixed state t > ^/a /ir is determined completely by the autocorre- 
lation function of random variable u(t). Randomsity of this variable in it's turn is connected 
to the chaotic motion of atoms inside of cavity. 
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